Structure and kinetics in the freezing of nearly hard spheres 
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We consider homogeneous crystallisation rates in confocal microscopy experiments on colloidal 
nearly hard spheres at the single particle level. These we compare with Brownian dynamics simu- 
ations by carefully modelling the softness in the interactions with a Yukawa potential, which takes 
account of the electrostatic charges present in the experimental system. Both structure and dy- 
namics of the colloidal fluid are very well matched between experiment and simulation, so we have 
confidence that the system simulated is close to that in the experiment. In the regimes we can 
access, we find reasonable agreement in crystallisation rates between experiment and simulations, 
noting that the larger system size in experiments enables the formation of critical nuclei and hence 
crystallisation at lower supersaturations than the simulations. We further examine the structure 
of the metastable fluid with a novel structural analysis, the topological cluster classification. We 
find that at densities where the hard sphere fluid becomes metastable, the dominant structure is a 
cluster of m = 10 particles with five-fold symmetry. At a particle level, we find three regimes for the 
crystallisation process: metastable fluid (dominated by m = 10 clusters), crystal and a transition 
region of frequent hopping between crystal-like environments and other (m / 10) structures. 



I. INTRODUCTION 

Crystallisation is a long-standing challenge, due not 
least to its local nature, where rare events on micro- 
scopic time- and length-scales initiate the phase tran- 
sition [1]. This lack of understanding of crystallisation 
can have very significant practical consequences, for ex- 
ample in control of drug production [2 . It appears chal- 
lenging to make much progress with conventional mate- 
rials, due to the local nature of nucleation events which 
lead to crystallisation, however particle-resolved studies 
of model systems such as colloidal dispersions which cap- 
ture the essential thermodynamics provide the necessary 
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FIG. 1: (color online) Radial distribution functions in exper- 
iment and simulation, (a) low density, (j) = 0.27, (b) higher 
density = 0.53. In both cases, simulation parameters were 
carefully adjusted to experiment. Line, simulation data, cir- 
cles, experimental data. 



detail required [3|. 

Colloidal 'hard' spheres are important in the under- 
standing of crystallisation. Few systems have received so 
much attention, not least because both simulations and 
experiments can access relevant timescales and particle- 
level structural lengthscales [4-8 . The general phe- 
nomenology of hard sphere crystallisation has been well 
established for a decade [9]: at low supersaturations, 
close to the hard sphere freezing transition at a volume 
fraction of (/)/ = 0.494 [10 , crystallisation is dominated 
by rare events leading to the formation of large nuclei. 
Higher supersaturation leads to a very strong rise in nu- 
cleation rate, and upon increasing the volume fraction, 
approaching the hard sphere glass transition, crystalli- 
sation has been observed at times less than the struc- 
tural relaxation time [|TT1 [12], while at higher volume 
fractions still (0 = 0.62), crystallisation is not seen on 
the experimental timescale. Despite this phenomenlogy, 
very large discrepancies have been found in nucleation 
rates predicted by simulation using biased ensemble av- 
eraging and experiment [H [5l [13^ , which remain unex- 
plained [m |15] . Neither the inclusion of polydispersity 
O [13] nor electrostatic charge [T6 in the simulations 
has resolved this situation, although the former has been 
shown to have profound and complex consequences for 
nucleation [T71[T8] . 

Charles Frank originally suggested that 13-membered 
icosahedra might suppress crystallisation in the Lennard- 
Jones system [T9|, and recently there has been a resur- 
gence of interest in the role of local structure in crys- 
tallisation. In simulations of hard spheres, five-fold sym- 
metry has been identifiedboth with the suppression of 
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FIG. 2: (color online) Clusters detected by the topological 
cluster classification. These structures are minimum energy 
clusters of the Morse potential with po = 25.0 [36 . 
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FIG. 3: (color online) Structural changes upon increasing den- 
sity in the nearly-hard sphere fluid. Lines are simulation, ac- 
cording to eq. [l] circles are experiment. Data for metastable 
fluids (some of which subsequently crystallise) are taken at 
times <^ Tx . Dashed lines are estimated freezing and melting 
volume fractions for our system, as described in section (ITB] 
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crystallisation [20l [21] and found at the centre of crystal 
nuclei [22]. Locally dense amorphous crystal precur- 
sors have been identified in the metastable hard sphere 
fluid [23] and have also been found in softened sys- 
tems [24l [25]. One of us identified a mechanism for 
crystallisation through increased crystal-like ordering in 
the fluid prior to the formation of a nucleus, thereby low- 
ering the free energy barrier [26-28 , and that this entails 
no change in local density [29]. It was also shown that 
in weakly size-asymmetric binary hard sphere systems, 
crystallites can form quickly, but apparently become 'poi- 
soned' lag. 

Pioneering particle-resolved experiments ^ identified 
local structure, and more recent experiments on 'hard' 
spheres too polydisperse to crystallise have shown a de- 
gree of fivefold symmetry which, along with local crys- 
talline order, has been related to slow dynamics [31]. 
Here we consider local structure in crystallisation in a 
particle-resolved colloidal model system. While such ex- 
perimental studies can in principle resolve mechanisms 
of crystallisation, quantitative comparison to simulation 
and theory is very challenging, due to the limited accu- 
racy with which colloidal volume fractions can be mea- 
sured [32], combined with the lack of control over (and 
often knowledge of) interparticle interactions upon which 
crystallisation rates critically depend |T6[ |33l [34] . Quan- 
titative agreement between experiment and simulation 
has been obtained in the case of heterogeneous crystalli- 
sation of nearly hard spheres, initiated by a wall, where 
the crystallisation rate is less sensitive to the volume frac- 
tion compared to homogenous crystallisation [35] . 

Here we present a careful comparison of experiment 
and simulation in a system of nearly hard spheres which 
undergo homogenous nucleation. We interpret our re- 
sults with a novel structural method, the topological clus- 



ter classification (TCC) [37H39] . which directly identifies 
a number of local structures. Our mapping between ex- 
periment and simulation reveals good agreement in crys- 
tallisation rates at the range of supersaturation we ac- 
cessed. We find that the metastable fluid is dominated 
by 10-membered fivefold symmetric structures reminis- 
cent of the 13-membered icosahedra proposed long ago as 
a mechanism for the suppression of crystallisation [19]. 



II. METHODS 
A. Experimental 

We used polymethyl methacrylate (PMMA) particles 
of diameter 2.00 /im with a polydispersity of 4.0 % as de- 
termined by static light scattering. This degree of poly- 
dispersity is insufficient to have much impact on phase 
behaviour [41 . Swelling of the colloids cannot be ruled 
out in the density and refractive index matching solvent 
mixture of cis-decalin and cyclohexyl bromide used [32] . 
However, as far as crystallisation is concerned, electro- 
static charge, which is not entirely screened by the tetra- 
butyl ammonium salt added contributes a further degree 
of uncertainty in determining the effective colloid volume 
fraction. If ignored, the effects of electrostatics are quite 
sufficient to leave measures of crystallisation rates quan- 
titatively meaningless [161 [M] • For this reason we map 
simulations carefully to the experiments. We use confocal 
microscopy (Leica SP5) to track the particle coordinates. 
Heterogenous nucleation is prevented by weakly sintering 
larger (3.5 /im) polydisperse colloids onto the wall of the 
sample cell. We imaged at least 50 /im from the wall and 
saw no sign of heterogeneous crystallisation. 
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B. Mapping simulation to experiment 

Crystallisation experiments are compared with stan- 
dard Brownian dynamics (BD) simulations, with a sys- 
tem size of N = 2048 and 10976 particles and a timestep 
of 0.1 simulation time units. Fitting of the experimental 
radial distribution function is carried out using Monte 
Carlo (MC) simulations with N = 2048 particles. Both 
the BD and MC simulations are carried out in the canon- 
ical (NVT) ensemble. Our particle-resolved experiments 
enable an innovation: simulations take as their starting 
point experimental coordinates sampled from a fluid at a 
time small compared to that required for crystallisation. 
These are treated with 160 MC sweeps to remove small 
overlaps resulting from coordinate tracking errors. Par- 
ticle interactions are modelled with a truncated Morse 
potential with a Yukawa component, which approximate 
the hard core and electrostatic charging of the colloidal 
particles respectively. 

I3u{r) = Petm [l + e^°(^-^/^^^) (^e^^^^-^/^^^') - 2^ 

^Pey (1) 

Here the shifted Morse component (left term) is trun- 
cated at r = (Jij (where it vanishes) and the Yukawa com- 
ponent (right term) is truncated at r = 2cr^j. {3 = l/ksT 
the thermal energy, cr^j is the mean of the diameters of 
colloids i and j and r is the center separation. The trun- 
cated Morse potential is fixed with strength em = 1.0 
and range parameter po = 25.0. The contact potential 
of the Yukawa contribution /Sey = Z'^ Xb /[{I -\- i^cr /2)'^ a]. 
Here Z is the number of charges on the colloid and is 
the Bjerrum length. The inverse Debye screening length 
is denoted hy hz = ^/^'n^Bpion where pion is the number 
density of (monovalent) ions. 

We fix the Yukawa parameters to the expermental 
data. Our approach follows [42, 43 , where the Yukawa 
interaction parameters jSey and ncj are adjusted such 
that the experimental radial distribution function is well 
reproduced by the simulation (Fig. [T]). In this case we 
find ncj = 30.0 ± 5.0 and jSey = 1.0 ± 0.25, which corre- 
sponds to a Debye length of 67 nm (or an ionic strength 
of 1.4 jiM) and colloid charge of Z = 200. These are 
comparable to previous work on similar systems [42H46). 

We treat polydispersity with a Gaussian distribution 
in (J with 4% standard deviation (the same value as the 
size polydispersity in the experimental system). The 
radial distribution function of each experimental state 
point was fitted for a (metastable) fluid with MC sim- 
ulation. We then quote the state point in units of 
^ = ypart/yhox whcrc the volume of the particles is taken 
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There are a variety of ways to estimate the freez- 
ing transition for a system of weakly repulsive spheres. 
Among the more accurate |34j appears to be to in- 
terpolate exact simulation results for hard-core Yukawa 



systems [40] with the hard sphere values. This yields 
volume fractions for freezing and melting for a weakly 
charged system. Since we use a slightly softened core 
here, we estimate the impact of this softening on the 
phase behaviour by calculating the Barker-Henderson ef- 
fective hard sphere diameter cTg// 



^eff 



j dr[l 



exp(-/3M(r))] , 



(2) 



where u{r) is the interaction potential, i.e. Eq. ([T]) 
or a hard core with the same Yukawa term. The ef- 
fective hard sphere diameters are <j^j^^^ = 1.021(7 
and crj^j^^^ = 1.018(7 for the hard-core Yukawa and 
the truncated Morse system we use here. To approxi- 
mately include the slight effect of the core softening, we 
scale the volume fractions for freezing and melting by 
{^efT^^ l^ln^^^f- The core softening then leads to 
a change of around 0.004 in (j) in addition to the effect 
of the Yukawa repulsion. Thus for our system, we esti- 
mate the freezing volume fraction ^/ = 0.487 and melting 
(j)rn = 0.537. Note that we neglect the effect of the 4% 
polydispersity, whose impact on the freezing and melting 
volume fractions we expect to be slight [41^ . 



C. Structural Analysis - Topological Cluster 
Classification 

To analyse the structure, we identify the bond network 
using a maximum bond length of 1.4(7 and a modified 
Voronoi construction [37^ . For bond lengths greater than 
1.4(7, the network in condensed systems is insensitive to 
the bond length. Having identified the bond network, 
we use the Topological Cluster Classification (TCC) to 
determine the nature of the local environment of each 
particle [37l |39] . This analysis identifies all the shortest 
path three, four and five membered rings in the bond 
network. We use the TCC to find structures topologically 
identical to clusters which are global energy minima of 
the Morse potential for the range we consider po = 25.0, 
as listed in [36] and illustrated in Fig. [5] 

Now the system we consider interacts not via a full 
Morse potential, rather our truncation takes the repul- 
sive component only, in a similar spirit to the approach 
Weeks, Chandler and Andersen (WCA) used for the 
Lennard- Jones model [47] . While that approach is well- 
known to reproduce accurately the fluid structure at the 
pair level, one might expect deviations for higher-order 
structure such as that probed by the TCC. In fact we 
found that for short-ranged systems, clustering is en- 
hanced in the case of truncation Unlike many 
analyses, for example those which use bond-orient ational 
order parameters [49], our emphasis on bond topology 
distinguishes between icosahedra and the 13- membered 
D^y^ structure illustrated in Fig. [2] which is the minimum 
energy cluster for the Morse potential with (po = 25.0). 
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We have also checked for the icosahedron and found smah 
quantities (^1%) [48 . In addition we identify the thir- 
teen particle structures which correspond to FCC and 
HCP in terms of a central particle and its twelve nearest 
neighbours. For more details see [37, 39 . If a particle is 
a member of more than one cluster, we take it to reside 
in the larger cluster. 



III. RESULTS 

A. Fluid structure: matching state point 

We begin our presentation of results by considering the 
stable and metastable fluids. For metastable fluids, the 
cluster populations are sampled from the period prior to 
crystallisation (see next section). The topological clus- 
ter classification (TCC) analysis shows a considerable in- 
crease in cluster populations as a function of volume frac- 
tion, as shown in Fig. [3] Some smaller clusters are sub- 
sumed into larger clusters for (j) > 0.55. We see a sharp 
rise in the fivefold symmetric 10-membered C2v cluster 
we term 'lOB' following [36 , such that by > 0.54, it is 
the dominant cluster in the fluid. We estimate freezing 
in our nearly hard sphere system at ^ = 0.487 as noted in 
section [II C[ We note that the experimental and simula- 
tion data are very well-matched in Fig. [3j This gives us 
confidence that the fluid structure of the experiments is 
accurately reproduced in simulation, in other words that 
the state point is well matched and that we are therefore 
in a strong position to investigate any possible discrep- 
ancies between the experimental and simulated system. 



B. Matching timescales 

The timescales in the experiments and simulations are 
matched as follows. The dynamics slow upon increas- 
ing and we define a structural relaxation time from 
the intermediate scattering function (ISF). We deter- 
mine the intermediate scattering function (ISF) F(k, t) = 
{cos(k-(r(t+t')— r(t'))^)) where the wavevector k is taken 
close to the main peak in the static structure factor and 
r is the location of each particle at time t and the angle 
brackets denote a statistical average. An ISF is shown 
in Fig. |4] We fit the tail of the ISF with a stretched ex- 
ponential F(k,t) = C ex.-p[{—t/ra)^] to obtain the struc- 
tural relaxation time r^- We obtain C = 0.99, = 30.5 
s, /3 = 0.8134. Tq, is then plotted as a function of <j) in 
Fig. [5] The data are well described by a Vogel-Fulcher- 
Tammann (VFT) relation = tq exp[I)^/(^o — ^)] 
where tq is a relaxation time in the normal liquid, D 
is the 'fragility index', and (j)o ~ 0.616 is the ideal glass 
transition volume fraction [28^. We assume the scaling 
of the dynamics at this volume fraction is the same for 
experiment and simulation. We equate the experimental 
and simulation structural relaxation times at (j) = 0.43. 



Thus both experiment and simulation Tq, are taken from 
the VFT fit (Fig. |5|. 

C. Structure in crystallisation 

The process of crystallisation is shown in Fig. [6j Here 
^ = 0.54. Note that some small crystalline regions are 
present at t = 0. We have previously demonstrated that 
even stable fluids have populations of particles in crys- 
talline environments [48 , as shown in Fig. |3j this rises 
markedly around the freezing transition. These images 
suggest that the metastable fluid shows relatively little 
change for 27.9tc^ < t, but crystallisation subsequently 
occurs at 27.9rQ, < Tx < 316.9rQ,. Henceforth we define 
i~x as the time at which 40% of the particles are iden- 
tified in crystalline environments. Moderate changes to 
this threshold have no impact upon our conclusions. 

We now turn to a TCC analysis of the crystallisation 
process. Figure [7| shows the population in each TCC 
cluster as a function of time for experiment (a) and sim- 
ulation (b). In both experiment and simulation, we iden- 
tify three regimes. For t < IOtq, in the experiment and 
t < 4:0Ta in the simulation there is little change in clus- 
ter populations. At intermediate times approaching Tx 
(here Tx = 107 and 64rc^ for this experiment and simula- 
tion run respectively), we see a steady growth in parti- 
cles identified in crystalline environments (predominantly 
FCC) at the expense of particles in fluid environments. 
Most notable of the non-crystalline clusters is lOB, which 
drops continuously throughout this period. At the level 
of this description, then, crystallisation is interpreted as 
the conversion of lOB clusters into FCC environments. 
At times larger than r^^, there is a further decrease in non- 
crystalline clusters. However, note that, on the timescale 
of these experiments and simulations, a reasonable pop- 
ulation (a few percent) of non-crystalline clusters remain 
at all times. 



D. Comparision between experiment and 
simulation 

We now compare crystallisation times in simulation 
and experiment. Recall that nucleation of 'hard' spheres 
is found to exhibit strong deviations between experiment 
and simulation [4, 14 . We compare crystallisation times 
as shown in Fig. |8] We see a reasonable agreement for 
moderate values of (j) > 0.56, but at lower supersatura- 
tion (j) < 0.55 or — (j)M ^ 0.01, we find an emergent dis- 
crepancy between experiment and simulation. While no 
mapping between experiment and simulation is perfect 
[32 , 34 , our careful analysis of state point and timescale 
leads us to believe that this discrepancy is not accounted 
for by a shift of (j). 

Now neither the simulations we employ here, nor the 
confocal microscopy experiments access the regime of low 
supersaturation where the formation of a large nucleus is 
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a rare event. In simulation, biasing techniques can be 
employed, and while similar methodologies are in prin- 
ciple possible in experiment [50 the kind of precision 
required to determine nucleation rates quantitatively re- 
mains some way off. An important point then, is that 
unbiased simulation and confocal microscopy experiment 
access similar regimes of super saturation. However the 
experiment has a rather larger system size than does the 
simulation. The simulation box size is typically 2000 and 
10,000 cr^ for the = 2048 and 10976 system sizes re- 
spectively, while the experiments are confined in capillar- 
ies of size 250crx250crx2500cr - 1.6xl0^cr^. The imaging 
volume (50(7 x 50cr x 25cr ~ 6.3 x lO^a^) is rather smaller 
than the whole system and crystals can nucleate outside 
this region (Fig. [6|. The rate of crystal growth has re- 
cently been determined in a very similar system [35 , and 
the associated timescales are of order IOO-IOOOtq,, sug- 
gesting that for long times, crystals can spread through- 
out the sample, so the relevant volume is that of the en- 
tire system, rather than just the imaging volume. Thus, 
decreasing the supersaturation to (j) — (j)M ^ 0.005, we see 
that experiments continue to crystallise, but for simula- 
tions, the time for crystallisation moves outside the ac- 
cessible timescale. Note that, at higher supersaturation, 
nucleation rates increase strongly, so the crystallisation 
time Tx is somewhat independent of system size. 

While the argument presented above is physically at- 
tractive, we seek a more quantitative validation, and for 
this we turn to the results of Filion et al. . Clearly, 
for the system to crystallise, at least one nucleation event 
must occur. Taking the nucleation rate for the high- 
est volume fraction in ref. [14 , for a monodisperse 
system J = 1.4 x lO^^a^r^ where tb is the time to 
diffuse a diameter in the dilute limit and is equal to 
O.63rQ,(0 = 0.43). The corresponding volume fraction 
relative to melting is <p — (j)M = 0.0051. For our system, 
this corresponds to one nucleation event every 4.0 and 
21.6 To,{(l) = 0.43) for N = 10976 and N = 2048 respec- 
tively. In other words, in our simulation timescales, we 
expect crystallisation. 

However, the rate of nucleation is highly sensitive to 
polydispersity [3 [13], which is 4% here. Recall this is ex- 
pected to have little effect on the equilibrium phase dia- 
gram [41]. While to the best of our knowledge, no precise 
predictions for nucleation rates in polydisperse systems 
have been made for the regime we access {<p > </>m), it 
seems reasonable to expect a shift of the nucleation rate 
as a function of (j) by ^ 0.015 which has been found at 
lower supersaturation \5\ [13] . We note that the shape 
of the particle size distribution, rather than just its sec- 
ond moment, can be important [17l[T8], however here we 
shall assume that the effect of polydispersity is to effect a 
shift in (j) of 0.015 in the nucleation rate. In other words, 
we take the rate for (j) — <pM = —0.01 for a monodisperse 
system to apply for (j) — (pM = 0.005 for our system. The 
nucleation rate at ^ — = —0.01 is four orders of mag- 
nitude lower than that at — 0m = 0.005 jlj. At such 
low rates, we expect no nucleation in simulation, but the 



larger system size in experiment (conincidentally four or- 
ders of magnitude larger than the simulation) is sufficient 
for nucleation to occur on our timescales, as is consistent 
with the crystalisation that we see. 

If this analysis is correct, for a monodisperse system 
at (/) — (j)M = 0.005 we should expect a much higher 
nucleation rate, and crystallisation on the simulation 
timescale. To verify this point, we carried out some sim- 
ulations with a monodisperse system. These are shown 
in Figs, [sj^a and b), and indeed crystallise in the regime 
of interest. We thus conclude that, in the regime we ac- 
cess, the discrepancy between experiment and simulation 
is likely due to the larger system size in the case of the 
experiments. 



E. Particle-level crystallisation mechanism 

Our single-particle analysis gives us the ability to shed 
some light on the mechanism of crystallisation. We have 
noted that the metastable fluid is dominated by the five- 
fold symmetric lOB cluster. In Fig. |9j we show the 
history of a single particle, throughout a simulation for 
(/) = 0.55 and (j) = 0.58 . We see that each particle fluctu- 
ates and is identified in a number of different structures, 
including local crystalline environments. Three regimes 
emerge: the metastable fluid, dominated by the lOB clus- 
ter, the final crystal, and a transition regime between 
the two. Note that in the transition regime, although 
the particle is often found in amorphous structures, these 
are rarely lOB. This suggests that lOB-crystal transitions 
may be somehow suppressed. This is consistent with 
long-standing ideas that five-fold symmetry can suppress 
crystallisation [I9l and very recent experimental j3T] 
work which suggests frustration between five-fold sym- 
metry and local crystalline order. We would thus expect 
that lOB clusters are rather stable (see below). During 
growth, however, it is possible that a crystalline surface 
may disrupt the five-fold symmetry in the fluid, lead- 
ing to more rapid transformation between lOB and crys- 
talline structures. Our finding of a transition regime fur- 
ther supports findings that crystallisation preferentially 
takes place in a region of high crystal-like bond orienta- 
tional order [261(281 [29 . 

We close by considering the stability of the fivefold 
symmetric lOB cluster. This seems to dominate the 
metastable fluid at densities where crystallisation occurs. 



In Fig. 10 we show transition probabilities from the lOB 
cluster to various geometries. We see there is a tendency 
to remain in the lOB cluster. That is, the lOB clus- 
ter shows a higher degree of stability than other clus- 
ters. In other words, for the nearly hard sphere fluid, the 
lOB is a locally preferred structure, similar to icosahedra 
and related polyhedra in glass-forming systems [49l [5T1 - 
|55] . This is consistent with previous reports that five- fold 
symmetry is favoured in 'hard' spheres [20l [21] [31] . The 
fact that this ordering tendency to structures of five-fold 
character such as lOB appears to be enhanced at higher 
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(j) is intriguing, as at higher volume fractions stih, 'spin- 
odal' crystahisation takes place in a small fraction of Tq, 
as found previously [TTJ [12]. In such unstable fluids, 
however, we cannot measure the lifetime of lOB clusters: 
the stability of lOB is exceeded by the thermodynamic 
driving force of crystal nucleation. 

IV. CONCLUSION 

We have carefully matched simulation to experiment 
for a nearly hard sphere system. For the regimes in which 
we can access crystallisation, the kinetics are similar in 
both simulation and experiment with the exception that, 
at lower volume fraction, experiments crystallise faster 
than simulations. We believe these are associated with 
the onset of low nucleation rates as the supersaturation is 
decreased. Under these conditions, the larger experimen- 
tal system size means nucelation events occur on acces- 
sible timescales, enabling crystallisation to be observed 
in experiments but not in simulations. While accurate 
simulation data has been obtained for monodisperse sys- 
tems across a wide range of (/), the effect of polydisper- 
sity characteristic of colloidal experiments has only been 
considered for nucleation rates too low for confocal mi- 
croscopy experiments to access. In order to be confident 
that no discrepancy exists, predictions for polydisperse 
systems in the regime accessible to experiments such as 
ours would be helpful. 

Our topological cluster classification reveals the mech- 
anism of crystallisation. In particular, around the freez- 
ing transition, nearly hard sphere fluids become dom- 
inated by a five-fold symmetric ten-membered cluster 
which we term lOB. By considering particle histories, we 
find that transitions between this lOB cluster and crys- 
talline environments are suppressed. Instead, after some 
time in a metastable fluid state, with occasional excur- 



sions to a crystalline environment, which usually occur 
through an intermediate structure, the particle finds it- 
self in a transition state, presumeably due to the prox- 
imity of a crystalline region which stabilises local crys- 
talline environments. In the transition state, the particle 
spends large amounts of time in a crystalline environ- 
ment, and almost no time in a lOB cluster, instead it is 
found in other amorphous clusters. Eventually, the par- 
ticle spends all its time in a local crystalline environment 
and is said to be crystalline. The transition regime we 
observe may be related to the 'cloud' identified in the 
case of softened particles [24| . 

Finally, we emphasise that, since absolutely hard 
spheres are not found in nature [34], it is essential to 
take account of the inherent softness in any experimen- 
tal system. However, comparison with true hard spheres 
suggests that the main effect of the softness we have con- 
sidered is to shift the state point such that we must con- 
sider effective volume fractions. This is consistent with 
previous observations that mapping the effective packing 
fraction to hard spheres results in a practically identical 
fluid structure [48] . 
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FIG. 4: Intermediate scattering function for experimental 
data at cj) — 0.43. The wavevector is taken close to the first 
peak in the static structure factor. Grey line is a streched 
exponential fit (see text for details). 




FIG. 5: (color online) Structural relaxation time as a function 
of (j). Solid line is a Vogel-Fulcher-Tammann fit (see text for 
details). 




FIG. 6: (color online) Confocal microscopy images of crystallisation in nearly-hard spheres for (j) = 0.54. (a) 600 s (2.3 ra), (b) 
4500 s ( 17.4 TO), (c) 7200 s ( 27.9 to) (d) 81900 s ( 316.9 to) bar=10 /im. 
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FIG. 7: (color online) Topological cluster classification analysis of crystallisation, experimental (a) and simulation data (b). 
Here 6 = 0.56. 
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FIG. 8: (color online) Crystallisation times in terms of Ta{(j) — 0.43) (a) and Ta{(j)). Circles are experimental data, dark and 
light squares are simulation data for polydisperse systems of = 10976 and N = 2048 respectively. Unfilled square is for a 
monodisperse system with N — 10976. Dashed lines are melting estimated as described in section al B| Solid lines are to guide 
the eye. Error bars extending upwards are lower bounds for crystalisation times determined from experiments (light lines) and 
simulations (dark lines) which did not crystallise. 
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FIG. 9: (color online) History of a single particle, (a) cj) — 0.55, (b) cj) — 0.58. Shaded areas mark the different regimes of fluid, 
transition and crystal, as described in the text. Data are shown from Brownian dynamics simulations. 
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FIG. 10: (color online) Transitions from the lOB cluster, (a) = 0.55, (b) cj) — 0.57. These are the probability for a particle 
to be found in a cluster at time t + rt, given that it was in a lOB cluster at t. Here n = O.OSr^ for both (a) and (b). Data are 
shown from Brownian dynamics simulations. 



